Bernoulli potential in type-I and weak type-II superconductors: Surface dipole 
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The Budd-Vannimenus theorem is modified to apply to superconductors in the Meissner state. The 
obtained identity links the surface value of the electrostatic potential to the density of free energy 
at the surface which allows one to evaluate the electrostatic potential observed via the capacitive 
pickup without the explicit solution of the charge profile. 



I. INTRODUCTION 

The Hall voltage is commonly used to measure the con- 
centration of charge carriers in conductors. In supercon- 
ductors this method is not applied, because of a miss- 
ing theoretical support. In this paper we show that the 
Bernoulli potential, which is closely related to the Hall 
voltage, can be used to the same end in superconductors. 

Any ideal collision-less electric fluid should exhibit a 
finite Hall voltage, but superconductors seem to escape 
this theoretical conclusion. The zero Hall voltage was 
first reported by Kamerlingh Onnes and Hof 1 already in 
1914. Though later analyses showed that their samples 
were in the mixed state what obscures an interpretation 
of their experiment, the zero Hall voltage was confirmed 
anyway. 2 

From a theoretical point of view, it was clear that 
there has to be a voltage balancing the magnetic pressure 
which acts on electrons via the Lorentz force. The report 
of Kamerlingh Onnes and Hof thus stimulated various 
speculations about the missing Hall voltage, see the crit- 
ical review by Lewis. 3 From various concepts proposed 
we mention the so-called contact potential - a potential 
step at the interface - which was expected to cancel the 
Hall voltage. Such a potential step might exist only if 
a charge dipole is formed at the interface of a supercon- 
ductor. The surface dipole we discuss here is a similar 
concept. 

The explanation of the zero Hall voltage turned out 
to be very simple. By contacts one monitors the differ- 
ences in the Gibbs chemical potential, often called the 
electro-chemical potential. 4 The Gibbs potential is com- 
posed of three components: the electrostatic potential, 
the kinetic energy and the correlation energy. None of 
these components is constant in the presence of diamag- 
netic currents, but their sum is constant in equilibrium 
in agreement with the observed zero Hall voltage. 

To eliminate the kinetic and correlation energies, Hunt 
proposed to access the electrostatic potential with a 
contact-less method called the Kelvin capacitive pickup. 4 
The first capacitive measurements appeared soon and 
they successfully proved the existence of a non-zero elec- 
trostatic potential. 5,6 



Both experiments 5,6 were done at temperatures well 
below T c , where the electrostatic potential has a simple 
form resembling the Bernoulli law, 7-9 
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To avoid confusion we note that the electrostatic poten- 
tial in equilibrium superconductors is called the Bernoulli 
potential for brevity, even if its actual form does not co- 
incide with the Bernoulli law. 

None of the early experimental data were sufficiently 
accurate to allow for a discussion of possible corrections 
to the plain Bernoulli potential (1). Nevertheless, the 
authors 5,6 made some conclusions in this direction and 
we find it necessary to comment on them in more detail. 



A. Bok and Klein 

Bok and Klein 5 claimed that their data agree with 
the plain Bernoulli potential (1). This conclusion has 
to be taken with reservations, however, because they 
measured the electrostatic potential as a function of 
the magnetic field B at the surface. They evaluated 
the velocity v of the superconducting electrons (briefly 
called the condensate velocity) from the London condi- 
tion mv = eA and the exponential decay A — Aoe~ x / x °, 
using B = V x A. At low temperatures, the London 
penetration depth depends on the density n of pairable 
electrons, Ag = m/(e fXon), therefore their experimental 
result can be expressed in terms of the magnetic pressure, 
enip = — £? 2 /(2/i ). 

According to the above arguments, the relation of the 
electrostatic potential to the magnetic pressure seems to 
be a consequence of the plain Bernoulli potential and the 
London theory. As we show below, the relation to the 
magnetic pressure is very general and holds also under 
conditions when neither the plain Bernoulli potential nor 
the London theory applies. 

In fact, Bok and Klein measured on indium which, be- 
ing a type-I superconductor, is not fully covered by the 
London theory. Moreover, they swept the magnetic field 
from zero to the critical value, while the plain Bernoulli 
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potential and the simple form of the London theory used 
above are restricted to low magnetic fields. At high fields, 
the condensate density at the surface is suppressed what 
results in v 4 and higher-order contributions to the elec- 
trostatic potential. 

Briefly, Bok and Klein have observed the magnetic 
pressure. But the link of their experiment to the 
Bernoulli potential (1) has to be taken with caution. 



B. Brown and Morris 

Brown and Morris have used a different setup which al- 
lowed them to achieve a much higher precision. 6 They did 
not control the magnetic field but the current in a thin 
wire. This current was scaled with the critical current. 
They announced in f 968 that their data reveal about 20% 
deviations from the screened Bernoulli potential, 
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discussed in more detail below. 

It should be noted that Brown and Morris expected de- 
viations which were predicted from the BCS theory in the 
same year. Adkins and Waldram had studied the elec- 
trostatic potential from changes of the BCS gap due to a 
current and they recovered the plain Bernoulli potential 
for zero temperature, while for finite temperatures they 
indicated a presence of additional contributions. 9 They 
were not, however, capable to derive these contributions 
in an explicit form or to estimate their amplitudes. 

Some corrections to the Bernoulli potential (I) were 
derived already before the BCS studies. Historically the 
first is the theory of Sorokin 10 from 1949 which covers the 
majority of effects recovered later. Although this paper 
is mentioned by London, 8 later it became forgotten. In 
1964 van Vijfeijken and Staas 11 took into account that 
the electrostatic field acts on normal electrons and ar- 
rived at the so-called quasiparticle screening. 

The quasiparticle screening is represented by the frac- 
tion of superconducting electrons n s /n by which eq. (2) 
differs from the simple Bernoulli law (1). In the exper- 
iment of Brown and Morris it accounts for 6% of the 
observed potential. In spite of its small magnitude, the 
quasiparticle screening is important with respect to the 
concept of the magnetic pressure. At finite temperatures 
one has to take into account that the London penetra- 
tion depth also depends on the density of condensate, 
A 2 = to/ '(e 2 Hong) = A§n/n s . Combining the screened 
Bernoulli potential (2) with the condensate velocity v 
found from the London theory, one finds that the elec- 
trostatic potential is temperature independent, i.e., it is 
given by the magnetic pressure with the density n of 
pairable electrons. Despite this importance of the quasi- 
particle screening, it should be noted, however, that the 
quasiparticle screening is not responsible for the 20% de- 
viations announced by Brown and Morris and discussed 
above. 



Corrections capable to explain the observed potential 
were first obtained by Rickayzen, 12 who showed that the 
electrostatic potential includes a pairing contribution, 
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The pairing term dominates close to the critical temper- 
ature T c , because n s = n(l — t 4 ), with t — T/T c , while 
n n = n — n s = nt 4 . 

According to (3), from ip close to T c one may deduce 
the density dependence of T c , which would be very attrac- 
tive with respect to designing new materials. Indeed, this 
important material property is otherwise deducible only 
from measurements applying a hydrostatic pressure or 
adding impurities to crystals. Unlike the later methods, 
a measurement of the Hall voltage or, respectively, the 
Bernoulli potential does not affect the electronic bands, 
the phonon spectrum or the electron-phonon interaction, 
therefore it offers a uniquely clear information about the 
material. 

All expectations were chilled by the next paper of 
Morris and Brown. 13 They admitted that deviations an- 
nounced in the first paper were due an incorrect estimate 
of the critical current and presented new highly accurate 
data for a wide range of temperatures. The observed 
electrostatic potential is perfectly equal to the magnetic 
pressure and exhibits no pairing contribution. They re- 
ported that this behavior is common to both type-I and 
weak type-II superconductors and for the magnetic field 
up to the critical value. 



C. Surface dipole 

The disagreement between theory and experiment re- 
mained unexplained for a long time and the question of 
the charge transfer in superconductors was left aside till 
the discovery of the high-T" c materials. For these lay- 
ered materials it was predicted ' 15 that the supercon- 
ducting transition induces a charge transfer from CuC>2 
planes to charge reservoirs. This transfer caused merely 
by the pairing mechanism has been confirmed by bulk- 
oriented experiments like the positron annihilation, 16 the 
x-ray absorption spectroscopy, 17 and the nuclear mag- 
netic resonance. 18 

Apparently, there are two groups of contradictory ex- 
perimental results. The pairing contribution is absent in 
the surface potential but a charge transfer is observed at 
internal interfaces in the bulk. 

As it was indicated recently, 19 there is a charge transfer 
at the surface which is the interface of superconductor 
and vacuum. This transfer forms a surface dipole which 
causes a step ipg in the electrostatic potential. The value 
of the potential step has been evaluated from the Budd- 
Vannimcnus theorem 20 
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where f c \ is the electronic part of the free energy den- 
sity (it does not include the electrostatic and magnetic 
parts). Rickayzen has obtained formula (3) from the gen- 
eral stability condition eip = — ^jn 1 and the free energy 
/ei = n s ^mv 2 . Using the same free energy one finds that 
the potential at the surface, <^o = + V>5, equals the 
screened Bernoulli potential (2). The surface dipole ips 
thus explains the observed magnetic pressure. 

In spite of the agreement between observed and the- 
oretically derived voltage one should be reserved about 
claims that the theory correctly describes the profile of 
the electrostatic potential in superconductors. Argu- 
ments against the theory are similar to those already 
raised in relation to the interpretation of the measure- 
ment of Bok and Klein. 

First, the measurement of Morris and Brown explores 
the entire range of magnetic fields from low up to critical 
values. The free energy employed by Rickayzen, however, 
applies in the limit of low magnetic fields only. 

Second, materials studied by Morris and Brown are 
type-I and weak type-II materials so that their behavior 
is not fully covered by the London theory. As we have 
shown recently, 21 even for low magnetic fields the electro- 
static potential depends on the Ginzburg-Landau (GL) 
parameter k. Rickayzen's formula is recovered for the 
extreme type-II superconductor k — > oo. For measured 
materials with n w l/\/2 the potential <p(0) at the sur- 
face is reduced by a factor 1/3 compared to Rickayzen's 
formula. 

Third, the surface dipole ips derived from the free en- 
ergy that covers only the low-field perturbation in the 
London approximation has the same shortcomings as 
Rickayzen's formula. Apparently, <p(Q) obtained from 
Rickayzen's formula can be far from the correct value 
and the same applies to ip$. Since the sum agrees with 
the experimental result, one can see that eventual errors 
tend to compensate each other in the resulting electro- 
static potential. In this sense, the surface dipole ip$ is 
consistent with the internal potential ip, since both are 
evaluated using the same free energy. 

D. Plan of the paper 

As demonstrated for Rickayzen's theory, the inter- 
nal electrostatic potential ip and the surface dipole ipg 
ought to be derived from the same free energy. We have 
shown 21 that for type-I and weak type-II superconduc- 
tors, the GL theory yields the internal electrostatic po- 
tential which is quite different from the one predicted by 
Rickayzen's formula. In this paper we derive the surface 
dipole within the GL theory. 

As in Ref. 19 we use the Budd-Vannimenus theo- 
rem. Here we employ this identity within Bardcen's 
extension 22 of the GL theory. Bardeen's extension offers 
two advantages. First, it naturally interpolates between 
the GL theory close to T c and the London theory at low 



temperatures. Second, it uses material parameters of the 
Gorter-Casimir two-fluid model, which have a transpar- 
ent density dependence. In contrast, the parameters of 
the original GL theory arc introduced in the limit T —» T c 
and T is replaced by T c wherever possible. Since T is an 
independent thermodynamic variable while j^j- =/= 0, one 
has to be careful when taking density derivatives. To 
evaluate the density dependence of the GL parameters, 
one has to recall the microscopic theory of Gorkov and 
take the density derivatives of the corresponding param- 
eters before the limit T — * T c is applied. 

The plan of the paper is as follows. In Sec. II we 
assume temperatures close to T c and an infinitesimally 
weak magnetic field. In this limit we derive a modifica- 
tion of the Budd-Vannimenus theorem which takes into 
account a non-zero charge density near the surface. In 
Sec. Ill we discuss a general system in the Meissner state 
and show that the formula for the surface dipole derived 
in Sec. II applies to any temperature and magnetic field 
below critical values. Section IV includes a summary. 

II. SURFACE DIPOLE 

Let us first estimate the thickness of the surface dipole 
from thermodynamic considerations. The pairing corre- 
lation is weaker on the surface than in the bulk, what 
results in forces pulling the Cooper pairs inside. Such 
forces are always balanced by the electrostatic field. The 
full understanding of this effect will require microscopic 
studies which are not yet feasible. From the BCS stud- 
ies it is known, however, that close to the surface on the 
scale of the BCS coherence length £ , the gap profile dif- 
fers from the value given by the GL theory. 23 We thus 
expect that the surface dipole is somehow linked to this 
'microscopic' modulation of the gap profile. 

To introduce the surface dipole on an intuitive level, 
let us assume that the system is close to the critical tem- 
perature. In this regime, the London penetration depth 
A and the GL coherence length £ are much larger than 
the BCS coherence length £o- Since the electrostatic po- 
tential induced by the diamagnetic current extends on 
scales of A and £ from the surface, 21 the surface dipole is 
very narrow on these scales. 

We can then define an intermediate scale L such that 
£o *C L <C £, A, as sketched in Fig. 1. On the scale 
L, the GL wave function changes only negligibly, i.e., 
ip(x) rts ip(x — > 0) = V(0) for < x < L. We note that 
the GL boundary condition d x tp = supports the slow 
change of ip close to the surface. 

As shown by de Gennes 23 , the GL wave function is lin- 
early proportional to the BCS gap, except for the surface 
region on the scale of the BCS coherence length £o- Fol- 
lowing microscopic theories giving the electrostatic po- 
tential in terms of the BCS gap, 14,15,24 we expect the 
electrostatic potential to have similar features, see Fig. 1. 

For £ <C x < L, the potential is well described by the 
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GL value <p(x). Since L <C A, £, the GL prediction of 
the electrostatic potential changes only negligibly in this 
region and it is convenient to introduce the extrapolated 
value tp(x) ~ tp(x — > 0) = <p(0). The extrapolation of 
the GL potential towards the surface, <p(0), has to be 
distinguished from the true surface potential ipo. The 
difference ips — tpo — tp(0) is caused by the surface dipole 
we aim to evaluate. 




FIG. 1. Electrostatic potential at the surface. On the scale 
of the London penetration depth A and the GL coherence 
length £, the electrostatic potential ip is linear in the GL wave 
function ip. The schematic short-dashed line in the expanded 
detail shows the deviation of the potential from the GL value 
on the scale of the BCS coherence length £o- The surface 
potential ipo is the experimentally observed value, while ip(0) 
is the internal GL potential extrapolated to the surface. In the 
insert, we show the virtual compression of the crystal lattice 
on the scale L, as it is employed in the Budd-Vannimenus 
theorem. The compression removes the lattice charge from 
the infinitesimal layer SL and correspondingly increases the 
charge density in the region framed by the long-dashed line. 



A. Budd-Vannimenus theorem at intermediate scale 

Close to the critical temperature, we can take the in- 
termediate region i ~ I as a homogeneous 'bulk' and 
follow the idea of Budd and Vannimenus. 20 ' 25 Let us as- 
sume a virtual compression of the crystal lattice such 
that the background or the lattice charge density is re- 
moved from the surface layer of an infinitesimal width 
SL. The perturbation of the lattice charge density in the 
infinitesimal layer < x < SL is Sp\ at = —put- The com- 
pression leads to an increase of the charge density in the 
layer SL < x < L, where 5 P \ a t = PutSL/L is selected to 
conserve the total charge. 

Now we recall the basic idea of the Budd-Vannimenus 



theorem. The lattice charge enters the jellium model of 
metals as an external parameter. If one changes this 
external parameter, the situation corresponds to doing 
work on the system, 



5W = S dx5pu t 



df 



' dp 
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S / dxSp\ at <p, 



(5) 



where / is the density of the free energy including the 
electrostatic interaction, and S is the sample area. Ac- 
cording to the Feynman-Hellmann theorem, the change 
of the electrostatic potential does not contribute to the 
work up to the first order in Sp\ at . Now we can proceed 
with the algebra. We split the integral into three parts, 
(0,SL), (SL,£o) and (£o,L). Since SL is an infinitesimal 
displacement, the potential in the layer < x < SL can 
be replaced by the surface value ipo. The surface region 
SL < x < £o gives a negligible contribution of the order 
of £o/L. In the remaining bulk region £o < x < L, the 
electrostatic potential is nearly constant and equals <p(0). 
The work thus reads 



SW = SSLp lat ( v (Q)- VQ ). 



(6) 



The work increases the free energy T of the system 



SW = ST = - 



d(USL) 
d(SL) 



SSL 
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SSL, (7) 



where f c \ is the spatial density of the electronic free en- 
ergy. Note that only the change of the electronic part 
in the 'bulk' is assumed. The surface energy does not 
change as the surface is merely shifted. The magnetic 
free energy changes negligibly, because the number of 
electrons in the layer L is not changed by the deforma- 
tion. Since the condensate velocity changes on the scale 
A ^> L, changes of the screening current vanish in the 
limit L/X — > 0. 

The first term in (7) results from the reduced volume, 
SL — ► S(L — SL), and the second one from the corre- 
sponding increase of the electron density, n — > n(l + 
SL/L). Equating (6) and (7) we obtain a modification of 
the Budd-Vannimenus theorem, 
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v(o)) = U-n%r- (8) 

on 



It describes the step of the potential at the surface due 
to the surface dipole in terms of the free energy. This 
relation is the main result of the present paper. In the 
remaining part we demonstrate how (8) can be used and 
show that it applies also at low temperatures, where the 
intermediate region cannot be defined. 

Formula (8) differs from the original Budd-Vannimenus 
(BV) theorem in four points. First, in the original BV 
theorem one evaluates the total potential step at the sur- 
face. Its values are of the order of Volts. Here we evaluate 
only the change of the potential step which appears as 
the system goes superconductive. The typical magnitude 
of this change is of the order of nano- Volts. Second, in 
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the BV theorem the surface potential is related to the 
potential ip^ deep in the bulk. In (8), the extrapola- 
tion of the internal potential towards the surface, <p(0), 
appears instead. Third, in order to cover systems at fi- 
nite temperatures, we use the free energy instead of the 
ground-state energy. Fourth, within the original Budd- 
Vannimenus approach, the electron charge density and 
the lattice charge density have to be equal because of 
the charge neutrality. In our approach, the density of 
electronic charge differs locally from the lattice charge 
density, en ^ —put, due to the charge transfer induced 
by the magnetic fields. 21 

In the limit of extreme type-II superconductors and 
weak magnetic fields one can compare formula (8) with 
the surface dipole evaluated in Ref. 19. In this limit, the 
free energy simplifies to / e i = — n s \mv 2 . The second 
term provides the current induced changes of the surface 
dipole derived in Ref. 19. Due to the bulk free energy 
foo one finds from (8) a finite potential step also in the 
absence of diamagnetic currents. 

For the purpose of the plot in Fig. 1, we have used a 
small magnitude of the potential step at the surface (the 
short-dashed line is rather close to the full line). In re- 
ality, even the current induced part of the step ip$ might 
achieve a magnitude much larger than the internal po- 
tential ip(0). A simple estimate of surface step for the 
experiment of Morris and Brown 13 ranges from few per- 
cent for T = 0.6 T c to about thirty times of the observed 
potential at T — 0.97 T c , see Ref. 19. Since we assume 
the limit T — > T c , the large values are more appropriate. 



/ = /ei + v(Piat + en) — ^-E 2 + — !— B 2 , (12) 

with the magnetic field B — V x A and the electric field 
E = -V<p. 

Variations of the free energy with respect to its inde- 
pendent variables A, ip,ip 7 n n yield the equations of mo- 
tion in Lagrange's form 



^ df df n 
-V — — + — = 0. 

dVv dv 



(13) 



For v = A the variational condition (13) yields the 
Ampere-Maxwell equation, for v = tp the Poisson equa- 
tion, for v = if) the GL equation, and for v = n n the 
condition of zero dissipation, 



eip 



dn 



(14) 



This condition allows one to evaluate the electrostatic 
potential in the bulk of the superconductor. 26 Of course, 
one can add any constant to the electrostatic potential. 

Formula (14) does not cover the surface dipole on the 
scale £cb therefore at the surface it provides the extrapo- 
lated bulk value <p(0). We can thus use (14) to rearrange 
the Budd-Vannimcnus theorem (8) as 



Piat^o = /el + <p(0) (plat + en) . 



(15) 



Now all terms on the right hand side are explicit quan- 
tities which one obtains within the GL theory extended 
by the electrostatic interaction. 26 



B. Application within the GL theory 

Now we demonstrate how the relation (8) can be used 
within the GL theory. To this end we introduce the GL 
free energy 

/ei = ^|HW-e*A)V| 2 + /cond, (9) 

where ip is the GL wave function, A is the vector po- 
tential, m* = 2m and e* = 2e are the mass and the 
charge of the Cooper pair, and / con d is the free energy of 
the condensate. It can be either the Gorter-Casimir free 
energy 

/cond = -\lTf~W \lT 2 ^l-%\ 2 (10) 

used by Bardeen 22 or the GL free energy 

/cond - C#| 2 + (11) 

The GL parameters a = j(T 2 -T 2 )/2n and = jT 2 /n 2 
depend on the temperature T and the electron density 
n = n n + 2|i/>| 2 , where n n is the density of normal elec- 
trons. Finally, we add the electromagnetic energy so that 
the free energy reads 



C. Convenient approximation 

In customary GL treatments, the electrostatic poten- 
tial and the corresponding charge transfer are omitted. 
For magnetic properties this approximation works very 
well, since the relative charge deviation, (put + en) /put, 
is typically of the order of 10~ 10 , leading to comparably 
small corrections in the GL equation. With the same ac- 
curacy one obtains the electronic free energy / e i. There- 
fore it is possible to evaluate the surface potential using 
the approximation 

<fo « -— (16) 

enoo 

which follows from (15) if terms proportional to 
(piat + en) / pi a t are neglected. By n^ we have denoted 
the asymptotic value of n deep in the bulk, i.e. the den- 
sity of pairable electrons, p\ at = —enoo. 

Within approximation (16) one does not have to eval- 
uate the potential profile and the related charge inside 
the superconductor. This is advantageous, in particular 
for systems of unknown material parameters dT c / dn and 
d^f/dn. 

Turning the argument around, from (16) one can see 
that the electrostatic potential at the surface cannot be 
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used to measure the material parameters dT c /dn and 
d"f/dn. This fact is already known from the experiment 
of Morris and Brown. 13 



III. MAGNETIC PRESSURE 

So far we have discussed systems close to the critical 
temperature, when the validity conditions of the GL the- 
ory are well satisfied. In many cases, however, the GL 
theory is used beyond the limits of its nominal applica- 
bility. In these cases the GL coherence length £ and/or 
the London penetration depth A are comparable to, or 
even shorter than the BCS coherence length £ so that 
the intermediate scale L cannot be introduced. 

It is possible, however, to follow the original formula- 
tion of Budd and Vannimenus and define L as the sam- 
ple thickness, i.e. L ^> £, A. In this case it is necessary 
to account for the energy of the magnetic field, since the 
infinitesimal compression 5L shifts the screening layer in- 
wards into the superconductor. 



A. Budd- Vannimenus theorem 

For L > the charge removed from the surface is 
placed deep in the bulk (the region of the scale of £, A 
gives a negligible contribution) so that the work on the 
charge reads 



5W = SSLiipoo - (^o)Plat- 



(17) 



Compared to the previous treatment we have merely re- 
placed the potential close to the surface by the value deep 
in the bulk. Similarly, the electronic part of the total free 
energy changes by 
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Finally, the shift of the screening layer by SL changes 
the magnetic energy by an amount STb = 
given by the magnetic pressure. Here Bq is the value of 
the magnetic field at the surface. 
From 6W = 5F e \ + 5T B follows 



Plat(</?0 - Woo) 



3. 
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As L is the thickness of the sample one can use the charge 
neutrality, p\ at — -en^. Since the value of the poten- 
tial deep in the sample is given by condition (14), using 
— gjj- from (19) one obtains 



etp, 



¥>o = - 



2p, en c 



/cl 



(20) 



The electrostatic potential observed at the surface is thus 
given by the magnetic pressure as observed by Morris and 
Brown. 13 



Note that deriving formula (20) we have not used 
many assumptions about the system. The condition of 
zero dissipation (14) is a general thermodynamic rela- 
tion. The Budd- Vannimenus relation (19), however, is 
limited to systems with a homogeneous jelly-like back- 
ground charge. This approximation is acceptable for 
conventional superconductors, where characteristic scales 
Cbcs, £ ; A are much larger than the elementary cell of the 
crystal. The applicability is questionable for the high-T c 
materials which due to the layered structure and a short 
coherence length are far from the jellium model. 

As noticed already by Bok and Klein, 5 there is a simple 
argument for the formula like (20). If one assumes a slab 
with magnetic fields Bl and Br on the left/right sides, 
the voltage difference gives 



Aat(<^0 



1 

2^o ' 



B R ). 



(21) 



The left hand side of this relation represents the electro- 
static force (per unit area) on the lattice, 



elst 



r-R 

J dx Ep lat = pi at (lf^ -ifo). 



(22) 



The right hand side is the Lorentz force i^or = BJ with 
the mean magnetic field B = \{Bl+Br) and the net cur- 
rent J — dx j given by Ampere's rule, B^—Br = poJ. 
Since the electrostatic field provides the only mechanism 
by which the force is passed from the electrons to the lat- 
tice, the two forces have to be equal, Fj^ OI = F e i st . This 
argument was, however, overlooked in the later studies. 



B. Test of the surface relation 

The Budd- Vannimenus theorem provides the electro- 
static potential (20) in terms of the magnetic field with no 
regards to the actual potential inside the superconductor. 
To link formula (20) with the more intuitive derivation 
from Sec. II, we show that (15) results in the surface 
potential (20) for any temperature. 

For the assumed geometry, the GL equation has an 
integral of motion, see Bardeen. 22 This integral can be 
obtained quite generally by the Legendre transformation 
of the free energy, 



= /-E 



(23) 



Indeed, if the fields v obey the equations of motion (13), 
the gradient V.g = vanishes, i.e. g = const. Deep in 
the bulk all gradients vanish, therefore g = /^°. 
From equations (9-12) one finds 



V ^-v v = — |v^ 



2 B 
Mo 



(24) 



G 



With the help of (23-24) and definition (12), one can 
express the electronic free energy as 

/e^/eT + ^IV^-^ + ^-^ + en). 

(25) 

At the surface, the GL boundary condition demands that 
= what implies E = 0. The free energy at the 
surface thus reads 

/ei = / e ? + ^ - P(°)(flat + en). (26) 

From (26) and the surface relation (15) it follows that 
Piat<Po = B 2 /(2[i ) + f°f. This value is identical to (20). 

Apparently, we can reverse the procedure. Starting 
from the general Budd-Vannimenus relation (20) and the 
general integral of motion (23), we can derive the surface 
relation (15). Accordingly, the surface relation holds for 
any temperature, provided that the free energy is a func- 
tional of the GL wave function and its first derivative 
only, / = f[tp, Vtp]. This functional can be an arbitrary 
one. 

Perhaps we should explain why we have derived the 
surface dipole from the Budd-Vannimenus theorem on 
the intermediate scale, although the more general deriva- 
tion from the integral of motion is available. There are 
two reasons. First, the intermediate scale provides at 
least a qualitative picture of the potential in the vicinity 
of the surface. This picture might be helpful if mea- 
surements sensitive to layers close to the surface will be 
designed. 

Second, within the intermediate scale the surface 
dipole is treated as a property of the superconducting 
condensate, what encourages us to hope that formula 
(15) or its approximation (16) can be used to obtain 
the surface potential also for cases when the magnetic 
field has a component perpendicular to the surface. In 
particular, we expect that it will be applicable to the su- 
perconductors in the mixed state, especially to evaluate 
the electric field generated by vortices penetrating the 
surface. 27 



IV. CONCLUSIONS 

In conclusion, the Budd-Vannimenus theorem was 
modified so that it is applicable to the surface of a su- 
perconductor. It allows one to evaluate the electrostatic 
potential on the surface from the free energy and the bulk 
electrostatic potential nearby. Formula (16) offers the ap- 
proximation of the surface potential from the free energy 
without the actual knowledge of the bulk potential. 

For plain surfaces we have recovered the experimen- 
tally established fact that the electrostatic potential 
equals the magnetic pressure divided by the density of 
pairable electrons. This experimental law was confirmed 



also for type-I and weak type-II superconductors, while 
the previous theoretical treatments were restricted to 
weak magnetic fields and extreme type-II superconduc- 
tors. The presented theory is free of these limitations. 

It was shown that thermodynamic corrections do not 
influence the surface electrostatic potential, measurable 
e.g. via contactless capacitive pickup. Consequently, 
contrary to earlier expectations, the density dependence 
of the critical temperature cannot be estimated in this 
way. On the other hand, the relation between the surface 
electrostatic potential and the magnetic pressure shows, 
that such a measurement allows one to determine the 
density of charge carriers without knowledge of any other 
material parameters. 

In this paper we have derived only the amplitude of the 
potential step. The detailed profile of the electrostatic 
potential including its modulation at the surface can be 
obtained by a microscopic approach like the Bogoljubov- 
de Gennes theory extended recently to cover the electro- 
static phenomena. 28-31 For microscopic calculations, the 
Budd-Vannimenus theorem can serve as a test of accu- 
racy of the numerical procedure, similarly as it is used in 
the theory of metal surfaces. 
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